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^ Abstract 

^ In a previous paper, we proposed a simple extension of the Minimal Super- 

h^- symmetric Standard Model which gives rise to thermal inflation and baryogenesis 

in a natural and remarkably consistent way. In this paper, we consider the A^ = 
^^ special case of our model, which is the minimal way to incorporate a Peccei-Quinn 

symmetry. The axino/flatino becomes the lightest supersymmetric particle with 
rria ~ 1 to 10 GeV and is typically over-produced during the flaton decay. In- 
terestingly though, the dark matter abundance is minimized for rria ~ 1 GeV, 
§^ /a ~ 10^^ to 10^2 GeV and \ijl\ ~ 400 GeV to 2TeV at an abundance coincident 

,^ with the observed abundance and with significant amounts of both axions and 

axinos. Futhermore, for these values the baryon abundance naturally matches the 

CN observed abundance. 
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1 Introduction 

In this introduction, we review some necessary background material on axions [1, 2, 3, 4] 
and thermal inflation [5, 6, 7, 8], and then give our motivation for the model presented 
in this paper. In Section 2, we present our model and describe some of its important 
features, including a detailed analysis of the flaton decay. Some technical calculations of 
flaton decay rates used in this section are given in Appendix A. In Section 3, we deter- 
mine the constraints on the model. In Section 4, we describe our numerical simulation of 
the leptogenesis and discuss the results. In Section 5, we summarize the physics of our 
model, compare it with related models, and discuss observational tests and signatures. 

1.1 Axions 

The main motivation for axions [9, 10, 11, 12] [13, 14, 15, 16] is to solve the strong CP 
problem [1, 3], though they are also one of the best candidates for dark matter [2, 4]. 
The strong CP problem arises because the QCD Lagrangian contains a term 

£. = »J^FF (1) 

which violates CP, and observations require 6 < 10^^'^ [17, 18] which is much smaller 
than would be expected naively. To solve the strong CP problem, Peccei and Quinn 
introduced a Peccei-Quinn (PQ) symmetry, U(1)pq, which is spontaneously broken by 
the PQ field 

where the pseudo-Nambu-Goldstone boson a is the axion. A U(1)pq-QCD anomaly adds 
an extra contribution to Eq. (1) giving 

Ce=(e- ^) -^FF (3) 

V faj 327r2 ^ ^ 

and, after the QCD phase transition, the axion relaxes to cancel 6, solving the strong 
CP problem. 

The axion decay constant fa is related to the PQ symmetry breaking scale by 

/. = ^ (4) 

A^ is the coefficient of the U(1)pq-QCD anomaly and is given by 

N= J2 P^ (5) 

«S{quarks} 

where the Pi are the PQ charges normalised such that (p has PQ charge one. The axion 
also couples to electromagnetism, indirectly via the U(1)pq-QCD anomaly and directly 
via the U(1)pq-QED anomaly. They generate a term in the QED Lagrangian 
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where C ~ 2 [19] and E is the coefficient of the U(1)pq-QED anomaly 

E = 2J2P^Q^ (7) 

i 

where the qi are the electromagnetic charges. The axion mass depends on the tempera- 
ture [2] 



0.1m,(0)(%^) for T»Aqcd 



where 



3.7 
D ^ 

maiT)c^{ - '\ T J — (8) 

ma{0) for T<Aqcd 

ma{0) ^ 6 X 10-^ eV f ^° ^''^ ] (9) 



fa J 



and Aqcd — 200 MeV is the scale of the QCD phase transition. 

Only invisible axion models with /„ much greater than the electroweak symmetry 
breaking scale are consistent with observations [20]. These models can be classified into 
two types: KSVZ [13, 14] and DFSZ [15, 16]. In the KSVZ model, the superpotential 
has a coupling 

W = \<j)XX (10) 

of the PQ field to new quarks x ^^^ X which become heavy after PQ symmetry 
breaking. In the DFSZ model, the superpotential has a coupling 

W = X,<P'H^Ha (11) 

between the PQ field and the Minimal Super symmetric Standard Model (MSSM) Higgs 
fields. 

Various astrophysical and cosmological observations constrain fa- Energy loss in 
Supernova 1987A gives a lower bound 



/,>109GeV (12) 

while axion cold dark matter abundance gives an upper bound. Axion cold dark matter 
is generated when the axion starts to oscillate coherently at the QCD phase transition 
due to an initial misalignment of the axion field, and also by the decay of PQ strings 
formed after PQ symmetry breaking. Assuming a randomized misalignment angle, as 
would be expected after a PQ phase transition, misahgnment generates [21] 

while the decay of PQ strings generates [21, 22, 23] 

r f \ '•''' 



compared with an observed cold dark matter abundance of [21] 



Q 



CDM 



0.2 



(15) 



However, if entropy is released after the QCD phase transition then the axions will be 
diluted ["' ■], see Section 3.6. 

The ADMX experiment [25] has searched for cold dark matter axions with axion 
mass in the range rria — (2 to 3) x 10^^ eV corresponding to fa ~ (2 to 3) x 10^^ GeV. 
For this range of axion masses, the resulting constraint on the ratio of the axion-photon 
coupling to axion mass translates to [19] 



to 0.8 < — < 3.6 to 4 



(16) 



1.2 Thermal inflation 



The main motivation for thermal inflation [5, 6] [24, 26, 27, 28, 29, 30] [31] is to solve 
the moduli problem [32, 33, 34], though it also solves the gravitino problem [35, 36] 
and may be able to explain the origin of the observed baryon asymmetry [7, 8] [37, 38] 
[39, 40, 41, 42]. Irrespective of these uses, thermal inflation is sufficiently natural that 
it might be expected to occur anyway. 

A flaton is a scalar fleld (p with negative mass squared at the origin and no quartic 
term 

V{(l)) = Vo-ml\(l)\' + ... (17) 

so that its vacuum expectation value 0o is much greater than its mass scale 



00 > m^ 
with m^ a typical soft supersymmetry breaking mass 

m,^ ~ ms ~ 10^ to 10^ GeV 



(18) 



(19) 



See Figure 1. 

Thermal inflation begins when the flaton is held at the origin by its flnite temperature 
potential and the energy density has dropped sufficiently for the potential energy at the 




Figure 1: Thermal inflation occurs when a flaton is held at the origin by its flnite 
temperature potential and Vq dominates the energy density. 




In a 



Figure 2: History of the universe with thermal inflation. The fractional density Qx = 
px/p of moduli, potential, flaton and radiation is plotted against the number of e-folds 
of expansion In a. 

origin to dominate. It lasts for about 10 e-folds and ends in a first order phase transition 
just before the temperature drops to the critical temperature 



Tc~ 



m^ 



(20) 



for the flaton to roll away from the origin. It is followed by a period of flaton matter dom- 
ination until the flaton decays leaving a radiation dominated universe at temperature 
Td. See Figure 2. 



1.2.1 Thermal inflation and the moduli problem 

Moduli are scalar fields with Planckian vacuum expectation values, and hence gravita- 
tional strength interactions. Their potential arises due to supersymmetry breaking, and, 
assuming supersymmetry breaking is transmitted to the observable sector via gravita- 
tional strength interactions, they would be expected to have vacuum masses of order 
the soft supersymmetry breaking scale in the observable sector 



''''inod 



m« 



(21) 



However, in the early universe, the finite energy density breaks supersymmetry. 
When H > rris this supersymmetry breaking dominates over the vacuum supersymmetry 
breaking and hence determines the moduh potential. When H drops below rris the 
moduh potential reduces to its vacuum form, but with the moduli typically displaced 
by a Planckian distance. The moduli then start oscillating with Planckian amplitude 
and immediately dominate the energy density of the universe, and, because of their 
relatively low mass and very weak interactions, they persist beyond nucleosynthesis 
with disastrous consequences [-32, -33, 34]. 

Infiation is typically invoked to rid the universe of unwanted relics, but one would 
expect an infiaton to have a mass > rris and hence primordial infiation to occur at a 
scale H > rris. However, the moduli are generated at if ~ rrig, and to a lesser extent 
by any phase transition at if < mg. Thus one wants infiation aX H ^ rris to dilute the 
moduli, but it is very difficult to realize primordial infiation at these scales. 

Thermal infiation [ ;, .>] on the other hand automatically occurs at if <^ rris- For a 
thermal inflation scale 

corresponding to 



V^'^ ~ 10^ to 10^ GeV (22) 



00 ~ 10^° to 10^2 GeV (23) 

thermal infiation provides enough dilution to rid the universe of moduli, but has a low 
enough scale not to regenerate them afterwards. Furthermore, it does not last long, 
only around 10 e-folds, and so does not destroy the primordial perturbations needed for 
structure formation, only shifting them to slightly larger scales. 

In more detail, following Ref. [6], the moduli produced before thermal infiation are 
diluted to an abundance 

and this is reduced further in the case of double thermal infiation, and those generated 
at the end of thermal infiation have an abundance 

^mod VoTd 



"^mod 



(25) 



with nucleosynthesis requiring [43] 



nmod < ^Q-12 (26) 



1.2.2 Thermal infiation and baryogenesis 

Thermal infiation provides a very natural solution to the moduli problem, but unfortu- 
nately is incompatible with most baryogenesis scenarios since it dilutes any baryons gen- 
erated before it and the temperature after fiaton decay is very low, typically (9(10 GeV). 
Fortunately, it gives rise to its own baryogenesis scenario [7, 8] [37, 38]. 
In Ref. [7], we proposed the simple extension of the MSSM 

W = XuQH^u + XdQHdd + KLHde + X^4>^HuHd + \k {LH^f + A^^XX + yA^^^ (27) 



with the key extra assumption that the soft supersymmetry breaking mass squared along 
the LHu flat direction is negative 



m 



LHu 



« - <J < (28) 



where m| and —m'jj are the soft supersymmetry breaking mass squareds of L and Hu- 
The flrst three terms in Eq. (27) are MSSM terms and the flfth gives neutrino masses. 
The term X^4>XX couples the flaton to the thermal bath. After thermal inflation, the 
X and X flelds acquire large masses 

M^ = A^0o (29) 

and so are not strongly constrained apart from that they should not damage gauge 
coupling unification. This coupling also induces |0| dependent renormalization group 
running of 0's effective soft supersymmetry breaking mass squared, tending to drive it 
negative at small |0|, as is required for a flaton. The term ^A^c/)^ stabilizes 0's potential 
at a value 

? (30) 

that turns X^cjP'HuHd into the MSSM /x-term with 

A, 






|/i| = A/,0c' ~ 



A„ 



m^ (31) 



The \^(f)^HuHd coupling also helps to ensure that the temperature after flaton decay T^ 
is high enough for dark matter to be generated. Thus we have a simple model of thermal 
inflation with the right scale for 0o if -^^ ~ ^pm and that also generates a /i-term of 
the right scale if A^ ~ A<^. 

Our key assumption, Eq. (28), seems at flrst sight dangerous since it implies a deep 
non-MSSM vacuum with LHu ~ (10^ GeV)^ and 

XdQLd + X^LLe = -fiLH^ (32) 

eliminating the /i-term contribution to LH^s mass squared [44]. To analyze the dy- 
namics induced by Eq. (28), we parameterize the potentially unstable flat directions 
as 

^=(000) , g = (^ ^/^^ U U ^ , d={d/V2 0) (34) 
= , x = , x = (35) 
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-10 
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Figure 3: Thermal inflation baryogenesis. 



The superpotential Eq. (27) reduces to 

W 
with the remaining D-term constraint 



11 11 

-Xdhdd^ + -Khde^ + X^(t)^Khd + -Xy {IKf + - 



D = |/lj2 _ |/,^|2 _ |;|2 ^ - |^|2 ^ - |g|2 ^ g 



^1 



The potential is 



V = Vo-m'M' + ml\l\' 



m 



,2|„|2 



2 I I, |2 I 2 I 7, |2 I 2| j|2 I 2| 

Hjhul +mHjhd\ +ma\d\ + mje 



(36) 



(37) 



+ I -AyXJ^hl + A^X^cP^Khd + -A^X^c^^ + -AdXdhdd^ + -A^XehdC^ + c.c. 

+ \Xulhl\ +\Xyfhu + X^(j?hd\ + \X^(f + 2X^(l)huhd\ + \Xdhdd\^ + \Xehde[ 



+ 



Xf,(t)^K + ^Xdd^ + -Aee^ 



+ i/[|/^„p-|/^,p-|/p + iMp + i|en (38) 



The dynamics is illustrated in Figure 3. Initially, all fields are held at the origin 
by their finite temperature potential. Once the energy density has dropped sufficiently, 
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Vo will dominate driving thermal inflation. After the temperature drops further, one of 
the unstable directions, (f) or Ihu-, will roll away from the origin. We assume Ihu rolls 
away first. Then ^AyXuf'h^ fixes the phase of Ihu, and \\ylh1\ and IX^t^hul stabilize 
its magnitude. The Ih^ field may partially reheat the thermal bath and so prolong the 
thermal inflation, but eventually will also roll away from the origin, ending thermal 
inflation. As rolls away, AuXu,4>'^huhd will force hd to become non-zero. \Xdhdd\ and 
\Xehde\ then constrain d and e to zero, shielding the dynamics from the dangerous 
non-MSSM vacuum in the direction of Eq. (32). Then A^X^cfy^huhd, and also ^A^X^cj)^ 
and the cross term from \X^(f)^ + 2X^(j)huhd\ , fix the phase of (p'^h^hd. As nears its 
minimum, the cross term from {X^l'^hu + X^(fP'hd\ rotates the phase of Ihu generating a 
lepton asymmetry, and at the same time \X^(j)'^hu\ gives an extra contribution to the 
mass squareds of Ihu and huhd, bringing them back in towards the origin. Thus we 
have a type of Affleck-Dine (AD) leptogenesis [15, 46, 47]. Preheating then damps the 
amplitude of the Ihu and huhd fields keeping them in the lepton preserving region near 
the origin [n]. The Ihu and huhd fields then decay, at a temperature in the MSSM sector 
above the electroweak scale, and their lepton number is converted to baryon number by 
sphalerons. Finally, the flaton decays, diluting the baryon density to the value required 
by observations, n^/s ~ 10~^°. 

1.2.3 Thermal inflation and axions 

The scale of the flaton vacuum expectation value, Eq. (23), coincides with that of the PQ 
fleld, Eqs. (4) and (12) to (15). This motivates uniflcation of thermal inflation and axions 
[5, 6], and several papers have studied this in detail [48, 49, 50, 51]. In order to introduce 
a U(l)pQ symmetry, Refs. [48, 49, 50, 51] extended the flaton sector, -^ 0i,02, to get 
a multi-fleld flaton axion model of the type studied in Refs. [52, 53, 54]. In this type 
of thermal inflation axion model, the main danger is that too many hot axions will 
be produced in the flaton decay to be consistent with Big Bang nucleosynthesis. This 
imposes signiflcant constraints on the parameter space. 

1.3 This paper 

The motivation for this paper starts from noting that our thermal inflation baryogenesis 
model, Eqs. (27) and (28), generates a very rich and complex but remarkably consistent 
cosmology from a very simple and constrained extension of the MSSM. 

As discussed in Section 1.2.2, the flrst three terms in Eq. (27) are pure MSSM terms. 
The fourth term generates the MSSM /i-term, forces Hd to become non-zero shielding 
the dynamics from the deeper non-MSSM vacuum, in combination with the flfth term 
generates the lepton asymmetry, couples the flaton to LHu and HuHd allowing the 
flaton to bring them back in towards the origin, and provides an efficient decay channel 
for the flaton allowing the temperature after flaton decay to be high enough for dark 
matter generation. The flfth term generates neutrino masses, stabilizes the LHu field 
during thermal infiation, and in combination with the fourth term generates the lepton 
asymmetry. The sixth term couples the fiaton to the thermal bath holding it at the 



origin during thermal inflation, and drives the flaton's mass squared negative at small 
101 . The seventh and last term stabilizes the fiaton field. 

Thus, all the terms we add, though simple and natural, do many different things. 
However, there is one exception. The term \X^(j)^ is simply introduced to stabilize the 
flaton's potential at the right scale. If we eliminate this term, the flaton's potential 
will in any case be stabilized by the |0| dependent renormalization of the flaton's mass 
squared induced by the coupling A^0XX- Furthermore, the model now has a U(1)pq 
symmetry and so automatically includes the axion without any need to introduce extra 
flelds. This simplest of axion models was studied long ago by Moxhay and Yamamoto 
[55] but has since been largely ignored. Its one disadvantage is that we have to tune 0o 
to the right scale, though this tuning is very modest. 

Thus we have an even simpler^ and more constrained extension of the MSSM, 
Eq. (27) with A,^ = 0, that generates an even richer and more complex, but as we 
will see, even more remarkably consistent cosmology. 

2 Model 

Our model ^ 



W = XuQH^u + XdQHj+ XeLHdC + X^^'^H^Ha + -X, {LH^f + A^0XX (39) 



1 
2' 

with the key parameter condition Eq. (28), is a simple extension of the MSSM incorpo- 
rating thermal inflation, baryogenesis and axions. In this section we describe some of 
its important properties. 

2.1 Fiaton potential 

The flaton potential has the form 

V{<P) = Vo - /(^«^ln^^±^'j mjl^r (40) 

where the function / encodes the |0| dependent renormalization of the flaton's mass 
squared. We set 

/(O) ^ 1 (41) 

to flx the deflnition of m^ at the origin and recover Eq. (17) 

V{<f)) = Vo-ml\<f)\^ + ... for 101 <m, (42) 

Away from the origin the potential simplifles to 

V{4>) = Vo -fL^ln^-^] m5|0|2 for |0| » m, (43) 



^It is arguable whether the absence of the 0^ term is more natural or not. 

^An obvious alternative would be to replace X^cji^H^Hd and A^ ~ -^pi^ with \^^4>^HuHd and A^ 



^GUT- 
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Now 



dV 






{2f + a^f')ml\<p\ 
(2/ + 3a^f + alf") 



Therefore, the potential has a minimuni at \(j)\ 

1 



We set 

to fix the definition of a^ at 
(46) and (47) give 



with 



(44) 
ml (45) 

/o = -2«0/o (46) 

/o = -1 (47) 

). Setting the vacuum energy to zero, Eqs. (43), 



Vq = -a^m^4>o 



and Eqs. (45), (46) and (47) give the physical mass squared at the minimum 



1 t/V 



m 



PQ 



2 t/|0p 



101=00 



"0/0 "0"i0 - «0"^0 



(48) 



(49) 



Note the suppression factor of a<^ in both Eqs. (48) and (49). In the vicinity of the 
minimum 

~1 



V{(l)) = Vo 



— a^ln 



^0 



+ ailn 



2 l„2 



^0 



2| J,|2 



^M 



(50) 



The value of 0o is determined by the |0| dependent renormalisation group running of 
the fiaton's soft supersymmetry breaking mass squared. The renormalization coefficient 
afj, is determined by the coupling \^4>XX i^i Eq. (39) 



«</, 



r-y$^|Axl'(^ + ^ + l^xl') 



■^n'^m 



(51) 



where m~^, m^ and A^ are the soft supersymmetry breaking masses and A parameter 
of X and x- Therefore, using Eq. (55), at least some of the Yukawa couplings should be 
unsuppressed 

El^xl'-l (52) 

X 

in order to obtain the correct scale for 0o- 

For definiteness, in our numerical simulation described in Section 4, we consider the 
simplest case 

fix) = l-x (53) 

We can then solve Eq. (46) to give 

00 = "^s exp ( ) (54) 



(55) 



To match Eqs. (19) and (23) then requires 

a^ ~ 0.05 
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2.2 Lightest supersymmetric particle 

In our model, the axino/flatino mass is generated by a XX loop containing the soft 
supersymmetry-breaking interaction A^A^^XX [55] 

"^'^^I^E^^x-ltolOGeV (56) 

X 

where we have used Eq. (52) to estimate A^,.. Therefore the axino will be the lightest 
supersymmetric particle (LSP) in our model. 

2.3 X ^iid X ^^d the PQ anomalies 

To maintain gauge coupling unification, x and x should form complete SU(5) multiplets, 
and so, assuming that they are not singlets, they should contain fields charged under 
SU(3). Thus the Peccei-Quinn field has couplings to both heavy quarks, contained 
in X a^^^d Xi ^^^d the Higgs fields H^ and Hd- Thus our model is a combination of the 
KSVZ and DFSZ axion models. 

Assuming that x can be decomposed into 1, 5 and 10 representations of SU(5), and 
X into corresponding 1, 5 and 10 representations, we can parameterize x and x as 

X = Nil + N55 + NiolO (57) 

X = Nil + N55 + N.olO (58) 

with the number of heavy quarks 

Ng = N, + 3Nio (59) 

Then Eq. (5) gives 

N = 6-Ng (60) 

with the existence of an axion requiring A^ 7^ 0, and Eq. (7) gives 

E = 12- ^N, (61) 

and therefore 

E 8 A 

N-3-N («^' 

As discussed above, x and x should form complete representations of SU(5) in order 
to preserve gauge coupling unification. However, although gauge coupling unification 
is preserved, the GUT scale gauge coupling becomes stronger than in the pure MSSM 
case. To preserve perturbative gauge coupling unification, we must restrict the size of 
the X and x representations. Using the parameterization of Eqs. (57), (58) and (59), 
this requires [56] 

N, < 6 (63) 
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2.4 Baryogenesis 

The potential derived from Eq. (39) is a simplified version of Eq. (38) 

V = V,-ml\<P\^ + ml\l\^-ml^\K\^ + mlJ^h,\^ 
+ { -A^XJ^hl + B\^<i?Khd + c.c, 

+ \9'{\K?-\h,?-\l\'Y (64) 

where, as described in Section 2.1, 




mim) = f\ -a^ln^^^-^ mi (65) 



and we have set d = e = since, as discussed in Section 1.2.2, they will be held at the 
origin throughout the dynamics. Despite setting A<^ = 0, the potential still contains all 
the terms needed for our baryogenesis scenario, and we expect the thermal inflation and 
baryogenesis to proceed as described in Section 1.2. The detailed differences will be due 
to the presence of the axion and the suppressed flaton vacuum mass, Eq. (49), both of 
which may affect the preheating and decay of the flaton. The PQ strings formed at the 
end of thermal inflation may also have some effect, possibly prolonging the leptogenesis. 
In Section 4 we describe and present the results of our numerical simulation per- 
formed in order to test in detail whether the baryogenesis proceeds as we expect, espe- 
cially that the preheating leads to a conserved lepton number, as we did for our original 
model of Ref. [7] in Ref. [S]. 

2.5 Flaton decay 

The flaton decays to Standard Model (SM) particles and axions 

r^ = TsM + Ta (66) 

2.5.1 Axion branching ratio 

The decay rate to (hot) axions is (see Appendix A) 



647r0^ 

The decay to SM particles is dominated by three channels: direct decay to SM Higgs, 
decay to bottom quarks via flaton-Higgs mixing, and decay to gluons via a XX loop 



SM — J- (t>^hh 



+ r^^feS + T^^gg (68) 
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The decay rates are calculated in Appendix A, giving 



'-'' - -/'"^-|-|'^VH^' j,_^^o{^,)+o{^i) 



r 
r 




%^ - O.la^AT^ (71) 

J- a 

where 7/, ~ 10~^ is defined in Eq. (247). For rrih ~ 125 GeV and rripQ ~ a^m"^, with 
a^p ~ 0.05 and m,^ ~ |/i| ~ mg ~ 10^ to 10^ GeV, we have roughly 

T^^hh r 10^ to lO'' for mpQ > 2m/, . . 

\ lO-^o 10-^ for mpQ<2mh ^ ' 

10^ to 10^ (73) 

10"^ to 10"^ (74) 

SO that either (p ^ hh or (p ^ bb dominate the SM channel with (p -^ gg always being 
negligible. Thus the SM axion branching ratio has the form 

^~i4=i^VmV/f:^| (75) 



Ta 




r^- 


*bb 


Ta 


r^^ 


99 



where 



and 



Ta \ '^A J V"^PQ/ V"^PQ 



/w^vr^+i^(i-f)' (76) 

e ^ ~ 0.02 (77 

ml 



see Figure 4. 

The branching ratio has a strong and somewhat complicated dependence on m\/mpQ, 
with the axion production suppressed if ^ 

^h ^ "^PQ - "<^"^J < ^l (78) 

The overall amplitude of axion production depends on the ratio 



f rnpq Y 



2 4 



which is highly suppressed in our model, see Eqs. (55) and (119). 



(79) 



•^This comes close to providing an anthropic explanation for the little hierarchy, since the Higgs 
thresholds could lead to sharp jumps in hot axion production, and hence via Big Bang nucleosynthesis 
in helium abundance, with anthropically relevant effects on stellar evolution. Unfortunately, the flaton- 
Higgs mixing channel makes this effect too small. 
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Figure 4: Fa/FsM versus mpg/m/j for m^ = 2|_B|, |/i| = 10^ GeV and m/i = 125 GeV. 



2.5.2 Flaton decay evolution 

While the flaton is decaying, its energy density evolves as [2] ^ 

p^ + "iEp^ = -T^ 

while the entropy increases as 

S _ FsmP</. 











S 


sT 


with 








3H^ 


= P^ + Pv 


and 








Pv = 


PSM + Pa 


Eq. (80) 


is easily solved 


to give 




P</,0( 


: a-^e-^^' 


while Eq. (81) gives 














s'/' = 


4 
3 


^27r2\ 
^45^ 


a' Jo 



I gl^^a^P^dt 
Jo 



(80) 

(81) 

(82) 
(83) 
(84) 

(85) 



*We neglect radiation from the AD sector. 
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Treating g^ as slowly varying gives 



PSM 



1- SM 



a^pA, dt 



and 



=4/3 



4 {2tx^ g,Y V 



45 

1/4 



SM 



1/4. 



a pd,dt 



(86) 



(87) 



Numerical solutions for p^, pr and gj (T^) /(7=„ (T) {S/ Sf) are shown in Figure 5. 

Flaton decay temperature There is no unique decay temperature, instead we define 
the flaton decay temperature T^ by 



or 



The time 



Pr(Td) = -Tl 



TT^ 1 

iTF^3*{Td)T^ = psM(^d) = 7^rsMr</ 



tn ^ 1.061 r 



-1 



corresponding to T^ is marked in Figure 5. At this time 

p^(td)^ 0.495 rj 
and after this time the entropy increases by a factor 

Si 



s^ 



1.946 



(88) 

(89) 
(90) 

(91) 
(92) 



Substituting the results of Section 2.5.1 into Eq. (89) gives 



Th 



87r4(7,(rd 
100 GeV 



mz 



\B\ 



m, 



l/^l 



1/2 , 

mpQ0o 



/I 



mt 



m 



PQ 



(93) 



,rv,2 I R|2 



m 



lOi^GeV 



U-M^^ 



J V 103 GeV y V "^PQ 



102GeV\2 



/I 



mt 



m 



PQ 



(94) 



which is plotted in Figure 6. 
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Figure 5: Flaton and radiation energy densities and entropy during the flaton decay: 
P(^/r|, Pr/r^, gj' {Ta) / gj (T) (5/ Sf) versus Ff^t. The dotted lines are the early time 

asymptotic solutions of Eqs. (100) and (101). The time corresponding to our definition 
of Td in Eq. (88) is marked in red. 
Td 




1 mh 



Figure 6: Flaton decay temperature versus fiaton mass: T^ versus mpq/mh for 0o 
10^1 GeV, |/x| = 10^ GeV, nih = 125 GeV and tua = 2\B\. 
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Early time decay evolution At early times, T^t <^ 1, 
therefore from Eq. (86) 



''^ ~ p^ oc a ^ 



5 



and 



Therefore, using Eq. (84), 



4 /27r2\ 4 



45 



^y^ (FsMi^)'^' oc a 



3/4 r. .-9/8 



(95) 
(96) 
(97) 



H 



PSM 




fld 



1.430 n 



Od 



. _3 . . _3 

r^(— ~o.69orJ — 

ad/ V^d, 



^y'(rd)^ ^^ s^fi^Y 



1/4 



^r'(T)^f 



TsMr,! 






Od 



0.828 TsMr,^ 



a 

Od 



0.750 



fld 



a 
ttd 
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(98) 
(99) 

(100) 
(101) 



3 Constraints 

Our model has a rich phenomenology generating many constraints. In this section we 
analyze the constraints on our model arising from the MSSM, neutrino masses, thermal 
inflation, baryogenesis, dark matter and nucleosynthesis. 



3.1 MSSM 

Stability of LHu in our vacuum requires 



l/ip > 2mlH^ 



(102) 



where frtLUu ^^ defined in Eq. (28). The MSSM /x parameter is generated by the fiaton's 
vacuum expectation value and has magnitude 



m 



I \ J,2| 



(103) 
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3.2 Neutrino masses 

The term \Xy{LHuY in Eq. (39) generates neutrino masses 

m^ = \KHl\ = \\^\v^sm^P (104) 

where v = 174 GeV is the electroweak symmetry breaking scale. The observed neutrino 
mixing [2i] 

Am^2 - 2.4 X 10-3 eV^ (105) 

Am^i ~ 8.0 X 10-^ eV^ (106) 



suggest neutrino masses 



and so 



ma ~ 5xlO"2eV (107) 

ma ~ 9 X 10-3 eV (108) 

mi < m2 (109) 



l^"l-«°«(iF^)to ("»' 



3.3 Thermal inflation and baryogenesis 

For the flaton and AD flat direction LHu to be unstable during thermal inflation, we 
require their soft supersymmetry breaking mass squareds to be negative 

-mj<0 (111) 

and 

- ^L„ < (112) 

The critical temperatures of the flaton and AD phase transitions are 

T, = ^ (113) 

and 

Tlh. = ^ (114) 



where 



pi=]yjK\" (115) 



X 

and 

PIh. = \{^\M' + ^92+91) (116) 

For the AD phase transition to occur before the flaton phase transition, we require 

n < Tlh^ (117) 
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therefore 

^ < ^ (118) 

PA PLH^ 



Combining this with Eq. (102) gives 



^0 2 , ^ai 

02 ^LHu < 2/52 

I LHu r^LHu 



< < l^rnl^< W^lf^l' (119) 



Now Eqs. (49), (52) and (55) give 



^PQ < ^ (120) 



and mpQ should be fairly close to this bound since we expect m^ ^ rUs ^ |/i|, while 
requiring m^ ~ mg > 100 GeV gives 

mpQ>20GeV (121) 

The roll away of the Affleck-Dine field LHu niay reheat the thermal bath and so 
extend the thermal inflation. This would dilute the moduli produced before thermal 
inflation by an extra factor [8] 

Aad ~ ^ (122) 

where we have taken the equality in Eq. (22) and hence Eq. (25) of Ref. [^], and noted 
that Fad *^ Hq so that the AD energy density is reduced by Hubble expansion rather 
than decay, and probably has a matter-like equation of state. 

After the AD phase transition, the AD field LHu settles to the minimum 

|/r = \hu\' = ll (123) 

of its potential 

\/ = yo+mi|/p-m^J/i„P+QA,A,/2/,2 + c.c.V|AXr+|A./Xr+^/ (|/i«r - If)' 

(124) 

'^ ^ m ^'''^ 



where ^ 



Eq. (110) gives 



For the flaton to be able to bring the AD field LHu back in to the origin at the end of 
thermal inflation, we require the flaton potential energy to dominate the AD potential 
energy 

a0"^j0o ^ ^IhJI (127) 



^The \LHu\ dependent renornialisation group running of to|^^ may reduce this value. 
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The extra dilution factor of Eq. (122) is sufficient to effectively remove the lower bound 
on 00 coming from requiring that thermal inflation sufficiently dilute the moduli, com- 
pared with the constraint 

00 > /o ~ 10*^ GeV (128) 

required by Eq. (127). The upper bound derived from Eq. (25) still applies, and using 
Eq. (48) gives 



VlO^GeVy 
The number of e-folds of thermal inflation is [6] 



<10i=^GeV /(^ ^'""d/A / m„,od y i^ lO^ GeV y |^ 10 GeV ^ 



(129) 



^-^inf zi^ I + Jlnf-^) -11 (130) 



6 \ in^.m^ I ' 3 



where the second term comes from Eq. (122). 

The baryon asymmetry generated in this model is [7, 8] 



r\-/ 



s V^ad/ \a^ml(f)l I \mLH^ 



V 10-2 J \ <po J \lGeVj V «W V ^u J \ m^ 



:i31) 



(132) 



Comparing with the observed value [21] 



!^ = 9 X 10-" (133) 



gives the constraint 



,12 r^ ^r 1 1 'f^L/n ad\ f ^d \ flO ^\ flO ^ eY\ frriLHu 



10'^ GeV ,/ ( ^^J ^^J ^— J [-^^ ^^J (134) 

Our simulation results. Figure 12 and Ref. [8], suggest that riL/nxD ^ 10^^ and so we 
get the rough upper bound 

00 < 10'' GeV (135) 

with values near this bound presumably being more natural. 
3.4 Axinos produced by the flaton decay 



The flaton dependent renormalization of the axino mass generates the effective radial 
flaton axino interaction 

rv-m- _ 

6ra + c.c. (136) 



(^aX^a r- ~2 

--= — or " 
V2 0O 
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where 

a-a ~ 10"^ (137) 

is calculated in Appendix A. 5. The flaton decay rate to axinos is 



Using Eq. (84), the flaton decay generates an axino number density 



(138) 



n-n 



9F C 9F 



and so, using Eqs. (89) to (92), a late time axino abundance 



a'p^dt = ^^^{e^^'-l)p^ (139) 



n 



a _ 2r^^aaSAe^^'^p^{U) _ 2.2 TdF^ 



(140) 



s mpql^ SfSd mpQTsM 

Current abundance Using Eqs. (89) and (138), the current axino abundance is 

n. - 5.6 X 10« (^) '^ (141) 

Therefore ^2^ < ^cdm — 0.2 requires 



m. < 1.8 GeV 



^V2 I 10 / llGeVy \ a-a \W^GeV 
which is at the lower end of its expected range Eq. (56). 



(143) 



3.5 Axinos produced by the decay of thermally generated NLSPs 

3.5.1 NLSP decay to axinos 

The superpotential coupling \^(j?HuHd leads to the decay of NLSPs to axinos via various 
channels with rate [58, 59] 

where we have neglected the decay products' masses, and A ~ 1 may contain a factor 
of m?z/m?^. 
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3.5.2 Axino abundance 

The thermal bath will generate NLSPs with number density 



nN = i, I SS^ (145) 



vr^ ./o /fc^ 



m 



exp^^^P^ + l 
/.(^) (146) 



whose decays will generate a late time axino abundance 

Ua 5d 1 / / a \ 



, . riNTN^adt (147) 



High Td limit In the high T^ limit 

Ua ^N-^d /"°° a^riN dt 



oo 



r- , / x^JN{x)dx 



-n-' Jo 



of 2 / \ T-^i/ 2 jy-i 2 



pl/2p 



rfgnmN)m% 



where we have used 



15C(5) 







Low Td limit Using Eq. (96), the axino production rate 



dlnT H T12 

and so, using 



V2^ V T J 

the maximum axino production occurs at 



21 
with our low Td limit corresponding to Td <^ T^r. 
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148) 



sirriN) Jo a{miqf 

T^%^N^d r f^^Y r (mN\dT 



(150) 
(151) 



x^ /7v(x) dx = ^^\ ' - 0.788 (152) 






^3/2^3/2 / m^\ 

rijv ~ — ^^ ^ exp ( — ^^] (154) 



T-Tn = ^rriN (155) 



We can analytically estimate the axino abundance in Eq. (147) using Eqs. (89), (92), 
(96) and (100), 



Sm)am?s(mH(TM)Jo \T J ^""Kt J T ^ ^ 



s 3S{T()a{TdYs{Ti)H{T, 

^SM 9* (^d) ^d ^N^a 

I '/ I ATI '/ I //'/ 

7r3(1.946) Jo 



Tl^^g^{TN)m% 



(157) 



4.4 X 10^ ^?^If| %^ f A ) (158) 



where we have used 

/ x^^ fN{x) dx = -{l- 2-") \Bi2\ vr^^ ~ 1.737 x 10 
Jo 3 



(159) 



Numerical A numerical solution of Eq. (147) can be obtained by writing it in the 
form 



r\^ 



.1/2 
SM 

where 



""'- 9f'iT,)gfHT^)^^^F,ix,) (160) 



V2 ^ „ o ^oo / X 3 



1 5 p^/^ ri 2 poo / \ J 

F,{x,) = g}{T^)g}{TN) ^^— / (-] nj,dt (161) 

versus 

^ 2^yum .^g2) 

'^ gj [Tn)Tn 

is plotted in Figure 7. 

For Xa -C 1, Eq. (158) gives 

Fa{xa)r^b.3xl (163) 

while for x^ :» 1, Eq. (151) gives 

F<j(x,) ~ 5.4 (164) 

Current abundance Using Eqs. (144) and (160), the current abundance is 

^^ - ^-^x^^liiv)" (^^^) 

" rf U''(^d)^*nT^);^10^GevAlGevJ»v 0o J'^'^^'^'^^ 

(166) 



24 



10° 



10 



10- 



-1 



10 



10" 



10- 



10- 



10- 



-3 



ji^a 












__ 




















^_ 






















_!_!_l^ 


!_!_L 








i::^:::- 


■ ■ 'Li: 


■ ■ ■ ■ 


































































































^ 


;L 


_ 


^ 












_ 


_ 


_ 
















.' 


^ 


= 


^X^ ^ 












^= 


^ 


^ 






























































































^ 




































y 









































^ 


= 














— 


^ 


^ 






























































































































/ 




































/ 































































— 


— 


— 


























































































































,y 










































































J 
























































































































































/ 






































/ 




































J 


/ 










^=^ 




















— 


— 


^= 




















































































































/ 






































/ 






































/ ^^ 












^ 


= 


= 














= 


= 


= 
















































































/ 






































/ 












_ 


_ 


_ 














_ 


_ 


_ 




^ 












= 


^ 


= 














^= 


^ 


^ 





























































































































































Xn 



10-1 



10° 



101 



Figure 7: Axino abundance versus flaton decay temperature: Fa versus Xa-, see 
Eqs. (160) to (164). F^, and hence the flaton decay temperature Td, is constrained 
by Eq. (167). 



Therefore ^2^ < i^coM — 0.2 requires 
F~Ax^ < 7.5 X 10-^ 



^rj/^ ( 9'J\T^)9'J\T^) \ ( lO'GeV \ /lGeV\ 
Aryj 1 103 \ mr, )[ m. 



J VlQiiGeV^ 

(167) 
wliich can be converted into a constraint on Xa, and hence the flaton decay temperature 
Td, using Figure 7. For xa -C 1, Eq. (163) gives 



(^.^0.19 4mViO^.*^/^(T.) 



t]/' \ 9fm 



N 



rriN 



WQeN) VlGeV/ 



'"^W12!1^V rH5ZiV (168) 



■ 



UlN 



and 



Td< 



m^f 



5'*(^J 



TV J 



25 L^*(Td)J 



ir^ /^<7y'(Td)^.'/'(T^)^ /^1 0^GeV \ /-IGeV 



Ari/2 

-•- SM 



103 



nia 



lOiiGeV 
(169) 
Note that the flaton decay temperature will be less than or similar to the neutralino 
freeze out temperature, which is about mi^i /20 [21]. Therefore we may require 



mpQ < 2mN 



(170) 
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to avoid direct production of neutralinos in the flaton decay. 

Axinos will also be produced by the decay of NLSPs after they freeze out. However, 
the standard Big Bang neutralino freeze out abundance is good match to the dark matter 
abundance, our freeze out NLSP abundance will typically be less than the standard 
abundance, and rria <^ wijv, therefore the axino abundance generated after freeze out 
should be safe. 

3.6 Cold axions 

3.6.1 Misalignment axions 

The axion mass turns on around T ~ Aqcd, see Eq. (8), and the axion starts to oscillate 
at a temperature To given by [2] 

m,(To) - 3Ho (171) 

with axion number density 

n„(To) oc m„(To) (172) 

and conserved axion number. The late time misalignment axion abundance is 

Ua 5* (To) na{Ta) 



Si s(To) 



(173) 



The abundance is suppressed by a factor A^ compared with the standard Big Bang 
abundance, where 

. _ {njsr _ Si m^T^^) / sJT^) _ g^T^) / To V'^ S^ 
^- njs 5(To) m,(To)/.(To) g..{T^^) \T^^ ) S{T,) ^ ^ 

Eqs. (8) and (171) give the standard Big Bang axion oscillation temperature 



rpU 
O 



,2- ' 



Si/BB 



(175) 



. 1.3GeV^l^)" (176) 

Low Td limit Our axion oscillation temperature can be analytically estimated using 
the results of Section 2.5.2 if the oscillation starts well before the flaton decay completes. 
For Td < To, using Eqs. (89) and (96), 

ttBB ('t^BB\3-7^4 
yj.7 ^ o { o ) o_ ^^^^^ 

Ho 
6^2 Tg^ g^ {To^)9* JTd) ^ bb\ 5.7 ^2 ^ yoN 

= wr?f — ^:m — ^ ° ^ ' ^''^ 
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and, using Eqs. (100) and (101), 



5(To) 



0.71 



r^*(^d)i 


4 


Psm(To) 

_PSM{Td)_ 



0.71 



g*{To) fTo 



9*iTd) \Td 



(179) 



Therefore, for T^ '^T^, the niisahgnment axion abundance is suppressed by a factor 



A, = 0.71 






SM 



g^T^) 



^,(TBB)^^(Td) 



3.7 . 3x4.7+1 

"'BB \ 7.7 



2x7.7 /J^ 



Td 



:i8o) 



High Td limit In the opposite limit, T^ ^ To, dilution due to the fiaton decay is 
negligible but the presence of hot axions produced in the fiaton decay increases H for a 
given T, and so reduces Tq, slightly enhancing the cold axion production. For Tj ^ Tq 



T 



and 



4.7 

fl-^^o) 



1 

2x5.7 



T 



BB 



9*iTl 



BB^ 



SM 

T7 



:i8i) 



:i82) 



Numerical We can numerically solve Eq. (171) using Eqs. (80), (81) and (82) to give 



r , \ 2X5.7 



r 



SM 






3.7 

^^" gl'\To)T^ 



Fo{Xa) 



where 



•^a 



' g*{Tf^) 
g*iTo) 



4x5.7 1/4 



Fo(x,)~1.0x' 



.TsM 
For Xq <^ 1, Eq. (178) gives 

while for Xa ^ 1, Eq. (181) gives 

Fo{Xa) ~ 1 

Eq. (174) then gives the suppression factor 



gj (Td)Td 
^y^(TBB)TBB 



0.26 



A. 



where Tf = min(Td, Tq) and 



6.7 



g*{To) 
g^T^"") 



4-7 1/4 

2x5.7 n ' 



g* [To) 



FaSxo) 



r , \ 2X5.7 



r 



SM 



g*{To) 



1/4 



gJ (To) To 
gl'\T^^)T^^ 



1/4 



gt'\To) Si 

g'J\T,)SiTo) 



(183) 

(184) 

(185) 
(186) 

(187) 



:i88) 
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Figure 8: Axion dilution factor versus flaton decay temperature: Fa^ versus Xa, see 
Eqs. (184), (187) and (188). 



is plotted in Figure 8. 

For Xa -C 1, Eq. (180) gives 



i"A.(xa)~0.69x, 



-1.96 



while for Xa ^ 1, Eq. (181) gives 



FaS^o) ~ 1 



3.6.2 String axions 

The axion abundance produced by the decay of PQ strings is [2] 

1 r° a^H"^ da 



Ua OC 



uj a 



(189) 
(190) 



(191) 



with the typical axion energy u (x H. Therefore, using Eq. (171), the late time abun- 
dance is 

n, 5(To) main) /"^° a'H da 



OC 



s S{ s{To) Jo alHo a 

Comparing with Eq. (174), our string axion abundance is suppressed compared with the 
standard Big Bang abundance by a factor 



^string 



"° a^H da 
«o-f^o a 



A misalign 



(193) 
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A string _ A misalign 

a r\ a 



string A misalign 



^strmg ^ ^ 



For Td < To, Eq. (95) gives 
while for T^ ^ Tq 

3.6.3 Current abundance 

Taking into account the dilution factor A^, Eqs. (13) and (14) become 

1.175 



a 



(194) 
(195) 



^a 



0.2 

a: 



Ja 



IQio to IQii GeV 



(196) 



For Td ^ T^^ ~ 1 GeV, the dilution factor A^ is slightly less than one, see Eq. (182), 
giving a slight enhancement of the cold axion abundance over the standard Big Bang 
case, which may provide an observational signal of our model in the future. For Td <^ 
T^^, Eq. (180) gives the dilution factor A^ ~ (Td/T^^^)^^-^^ and so, using Eqs. (4), 
(94) and (176), we have fi^ oc /-°-^^ 

3.7 Hot axions 

Hot axions will be produced in the flaton decay with energy density 



P, 



hot 



PSM 



4/3/ 



^.(Td)^.T(T) 



SM 



(197) 



where Fa/FsM is given by Eq. (75) and plotted in Figure 4. The axion production is 
highly suppressed in our model due to the suppressed flaton mass, but will typically be 
higher than the thermally produced axion energy density [2] 



P'^ 



10" 



10^1 GeV Y 



PSM V /a 

The hot axions produced in the flaton decay have current momentum 

a mpQ 



Pa 



flo 2 



(198) 



(199) 



where a is the scale factor at the time they were created and oq is the scale factor now. 
In particular, the current momentum of an axion produced at td is 



Pd 



OdrripQ 
flo 2 

Sy'gli\T,)T, 
1.48 X lO^^eV 



mpQ 
2 



1/3 



gj (Td)Td 



(200) 
(201) 

(202) 
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so, comparing with Eq. (9), they are marginally relativistic now. The current number 
density spectrum is 



Pa^ =(- ^^^ = °r Vr 3 (203) 

which may provide an observational test of our model in the future. Assuming that the 
hot axions are still relativistic now, their current energy density is 



fi!j°* ~ 2 X 10'^ 



100 

,9*iTd) 



3 F 

(204) 
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3.8 Nucleosynthesis 

The fiaton decay should complete before Big Bang nucleosynthesis, requiring 

Td > 10 MeV (205) 

The good agreement between the Standard Model prediction for the ^He mass frac- 
tion [uj] 

y/^ = 0.248 (206) 

and the observed value [21] 

Y°^' = 0.25 ± 0.01 (207) 

constrains any additional contribution to the energy density at the time of nucleosyn- 
thesis from hot axions produced in the fiaton decay [61] 



Pa 



PSM 

Eq. (197) then gives the constraint 

r 



< 0.14 (208) 

BBN 



- < 0.3 (209) 



J- SM 

which is automatically satisfied in our model, see Eqs. (75) and (120), and Figure 9. 

3.9 Summary of constraints 

Figure 9 shows the constraints on |/i| and ttipq for fa = 10^^ and 10^^ GeV, while 
Figure 10 shows the constraints on fa and |/i| for mpq = 0.1, 0.05 and 0.025|/x|. There 
are many strong constraints but nevertheless interesting regions of parameter space 
survive. The axino mass is required to be at the lower end of its expected range, 
rria ~ 1 GeV. Avoiding over-production of axinos requires /„ > 10^^ GeV, which nicely 
coincides with the more natural range of parameters for our baryogenesis scenario near 
the rough upper bound fa < 10^^ GeV. Avoiding over-production of axions splits the 
allowed parameter space into two regions. For fa ~ 10^^ GeV, the axion abundance is in 
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Figure 9: Parameter space constraints, |/i| versus mpq, for 00 = 3/a/\/2, mh = 125 GeV 
and niA = 2|-B|. Top fa = 10^^ GeV and bottom /„ = 10^^ Q^V. ■ hot axion over- 
abundance at Big Bang nucleosynthesis, Eq. (209). thermal inflation baryogenesis 
consistency constraint, Eq. (120). flaton mass constraint, Eq. (121). I axion cold dark 
matter over-abundance, Eq. (196) for misalignment or m high string estimate. 

■■ axino cold dark matter over- abundance: ■ Eq. (142) for «„ = 0.1 and solid m^ = 
1 GeV or line rria = 2 GeV; ■ Eqs. (167) and (170) for nia = 1 GeV and solid ttin = 
200 GeV or line niN = 100 GeV. 
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good agreement with the observed dark matter abundance, but to avoid over-production 
of axinos the MSSM /i parameter is confined to the range 700 GeV < |/i| < 2TeV. For 
fa ~ 10^2 GeV, the axino constraint weakens, but to avoid over-production of axions the 
decay temperature has to be low enough to dilute them, which requires lower values of 
|/i| in the range 400 GeV < |/i| < 800 GeV. 

The tightness of the dark matter over-production constraints implies that the dark 
matter should be composed of significant amounts of both axions and axinos. For 
example, Eqs. (142), (168) and (196) give 
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fi^ oc 0o^r; (7^d«7) (211) 


na oc 0^-^^^ (Td>lGeV) (212) 


therefore the dark matter minimum at fa ~ 10^^ GeV, /x ~ 1 TeV has 


nt = 7n^ (213) 


na = 1.7 Vta (214) 


and 


0.^/6 \°''^'' / a- \0-65 


r 102 1 


0.09 


r 102 1 


0.14 / -n \ 0.14 
1 't' \ 


L^*(^d)J 


_g*{TN)_ 


K'^smJ 


^^ / nia yo2 /200GeVy-^^ / A 


N 0.63 / 

Mevy V 


so/SHi \ 


' VI GeV/ \ niN ) V200 


5.6 X 108 GeV" V '"^"' 


Alternatively, Eqs. (132) and (196) give 


fib oc T^<pf (216) 


na oc 0j-52y|9^ (Td«lGeV) (217) 


therefore, fixing fib, the dark matter minimum at fa ~ 10^^ GeV, /i ~ 600 GeV has 


fia ^ fif = 1. 
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,5.6 X 108 GeV" 

To need to be at a minimum of the dark matter abundance may seem like tuning, 
but, taking into account anthropic selection effects, it can be quite natural for our local 
part of the universe to be found at such extrema of parameter space. Similarly, our 
mixed (supersymmetric!) dark matter may also seem odd, but is due to our extremum 
of parameter space, and follows the observed trend of an over-complicated composition 
of our local part of the universe, which is again well motivated anthropically. 
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4 Numerical simulation of the leptogenesis 

We closely follow the methods of Ref. [8]. We performed a three dimensional lattice 
simulation with periodic boundary conditions using the algorithm described in Ref. [8] . 
The potential was that of Eq. (64) with / given by Eq. (53), the D-term 

£) = |/i„P - |/irf|2 - |/|2 = g2 (220) 

and gauge 

constraints exactly conserved by the algorithm, and canonical gauge invariant kinetic 
terms. The numerical parameter e was introduced to cut off the singularity at hu = 0. 
The physical parameters and fields were rescaled by a typical soft supersymmetry break- 
ing mass m and a typical flaton expectation value Mpi or AD (/, hu or hd) expectation 
value Mad, as described in Refs. [7, ]. As the initial condition, we used the A = 4m 
case from Ref. [8]. 

4.1 Simulation parameters 

The lattice volume L^, number of lattice points A^^, time step At, and D-term constraint 
singularity cutoff e, were taken as ^ 

L = 200m-^ , A^ = 128 , At = 4xlO-^m-^ , e = 5 x lO^HfAD (222) 

We tested our results using different values of these numerical parameters. Limited 
computing power constrained us to A^ < 128. 

The fc-modes allowed by this lattice are k = \/k\ ^k\^k\ with 

ki = ^ (223) 

where rij = 0, 1, . . . , A^/2 with rij = 0, A^/2 having degeneracy one and the rest degener- 
acy two. This spans the range 

0.033 m= — <k< ^ — = 3.5 m (224) 

The physical parameters were taken as follows ^ 

m5 = 0.217m2 , a^ = 0.051758 , ml = ml (225) 

m^^ = 1.510m2 , m^^ = 3.533 m^ , m^ = 1.323 m^ (226) 

^It was necessary to take L larger than in Ref. [^] due to the suppressed flaton mass of Eq. (49) . 
'^The mass squareds were set to three decimal places to avoid any accidental resonances and also to 
fit the MSSM constraints described in Ref. [7]. 
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\A^\= 0.25m , |5| = 1.15m (227) 

|A^|=0.2mMAD , \X^,\ = 2.0mM^^ (228) 

Mad = 10~^Mti (229) 

arg {A,X,) = arg (BX^) = (230) 

Note that the phases of Eq. (230) are not physical and can be adjusted by field rotations. 
The CP phase was taken as 

(n-^ CP+ 

20 
arg(-5M^) = <^ n CPO (231) 

n+^ CP- 

\ 20 

For the choice of parameters above, we have [8] 

0o = O.7OMti , 0c = 0.6700 , /o = 1-5 x 10"^ 0o (232) 

where |0| = 0c is the critical value at which the niinimum of the AD sector switches to 
the origin 



'c \ _^H^- ^L 



4 

rn-, — rn.- 

,12 (233) 

The flaton mass squared eigenvalues at = 0o are 

m^Q = a0mj = 0.011 m^ , m^ = (234) 

and the Affleck-Dine sector mass squared eigenvalues at = 0o, / = /?-„ = /i^ = are 

mlH^=0.37m'^ , m^"^^ = 0.54m2 , m^+^jj^ = 3Am^ (235) 

and 

|/i|2 = 0.94m2 (236) 

Characteristic potential values are 

Vo = ^a^ml4 , Vi = 5.9 x 10"=^ Vq , ¥2 = 3.2 x 10^^ y^ (237) 

where 

V, = Vo- V{0, Ko, 0, h) (238) 

is the depth of the Affleck-Dine sector minimum when the fiaton is still at the origin, 
and 

V^2 = V^(0o,/i„o,O,/o) (239) 

is the height to which that point is lifted when the fiaton reaches its minimum. Here 
V^(0, hu, hd, I) is the potential. 
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Figure 11: Magnitude squareds of the flaton and AD fields averaged over the lattice as 
a function of time (CP—): (|0P)/0O' (I^P)/</'o' (K*rfl")/*/^o versus mt. 



4.2 Results 

The overall dynamics can be seen in Figure 11. Initially, the fiaton and the Higgs 
field hd are near the origin, the AD field / is near its thermal inflation minimum at 
\l\ = Iq, and hu is determined by / and h^ via the D-term constraint. As (p rolls away 
from the origin, ending thermal inflation, it first forces h^ to become non-zero, and then, 
as it crosses |0| = 0c, it causes the AD potential to flip, forcing / and hd in towards 
the origin. The AD fields then oscillate about the origin while the flaton rolls beyond 
its minimum and then back in towards the origin. Note that the flaton evolves slowly 
relative to the rapid AD field oscillations since the flaton mass scale mpq ~ a J rriff, is 
suppressed relative to the AD sector mass scales. As the flaton returns past |0| = 0c, 
the AD potential flips back to its original form, forcing / back out to |/| ~ Zq where it 
oscillates rapidly. The flaton nears the origin and then rolls out again, crossing |0| = 0c 
and flipping the AD potential for the final time, bringing the AD fields back in again 
towards the origin. The flaton then settles to oscillate about its minimum since the 
build up of gradient energy (preheating) has by this time drained enough kinetic energy 
from the homogeneous mode to prevent it from returning to |0| < 0c. The AD field 
/ initially oscillates around the origin with large amplitude but gradually settles down 
closer to the origin. 

The important outcome of this dynamics is shown in Figure 12. The AD field / 
is initially sitting in its thermal inflation minimum, and hence the lepton number is 
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Figure 12: Lepton number density averaged over the lattice as a function of time: 

mLHAnL)/V2 (CP+, CPO, CP-) versus mt. 

zero. The flaton then brings it in towards the origin with rotation, as described in 
Section 1.2.2, and hence lepton number is generated. This is initially conserved while 
/ is held in the symmetrical lepton number conserving part of its potential near the 
origin. Then the flaton forces / back out towards its lepton number violating thermal 
inflation minimum, about which it oscillates, and hence the lepton number oscillates 
violently. The flaton then again brings / back in towards the origin with rotation, and 
the corresponding lepton number is initially conserved as / is held near the origin by the 
overextended flaton field. As the flaton settles down to oscillate about its minimum, the 
/ fleld is held less tightly and spreads out to feel the lepton violating outer parts of its 
potential, and hence the lepton number decays. This decay is halted as / settles down 
closer to the origin, leaving a residual conserved lepton number. Although the dynamics 
is very complicated. Figure 12 shows that the lepton number is controlled by the CP 
violating phase in our potential, Eq. (231). 

The flaton dynamics is shown in more detail in Figs. 13 and 14. The flrst oscillation 
of the flaton is essentially homogeneous, but the angular dispersion becomes signiflcant 
in the second oscillation, limiting its amplitude. The radial flaton then settles down to 
oscillate about its minimum, with fairly small amplitude and dispersion, as can be seen 
in Figure 14. The axion on the other hand settles down to large amplitude oscillations 
in its dispersion, as can be seen from the late time oscillations in the mean squared and 
variance of the complex flaton fleld in Figure 13. 

The effect of the flaton dynamics on the AD sector is shown in more detail in Fig- 
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Figure 13: Mean squareds (solid) and variances (dotted) of the complex and radial flaton 
fields as a function of time (CP-): m?/<Pl (|0-(0)P)/0g, m)V<Pl ((l0|-(l0l))')/0^ 
versus mt. 
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Figure 14: Mean displacement from minimum squared (solid) and variance (dotted) of 
the radial flaton field as a function of time {CP—): (|0| — 0o)^/0O' ((I*/*! ~ (l</'l))^)/0o 
versus mt. 
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Figure 15: Magnitude squared of the flaton field and mean squareds (solid) and vari- 
ances (dotted) of the AD fields as a function of time {CP-): (|(/)p)/20g, |(/)|V^o 
(1^ - {l)?)/ll \{Ki)\Vll {\Ki - {h,)\')/ll versus mt. 



ure 15. The first two large amplitude fiaton oscillations force the AD sector in towards 
the origin, then out, then back in again. As the fiaton settles down, the AD fields 
mean squareds quickly settle to the origin but the variance now becomes significant. 
The variance of / is initially large but importantly decreases over time leading to the 
conservation of the lepton number seen in Figure 12. The variance of h^ on the other 
hand is initially small but grows over time. 

The preheating of the fiaton and AD sectors is shown in Figure 16. The initial roll 
out of the fiaton gives rise to tachyonic and angular preheating, exciting modes with 
k < m,/), as can be seen in the top left graph. Thereafter, the fiaton preheating becomes 
broader, but is still mostly limited to modes with k < m^, and the spectrum settles down 
to a fairly stable distribution well within the cutoff of our lattice, as can be seen in the 
bottom left graph. The AD preheating is more complex, as discussed in Ref. [8], and 
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Figure 16: Kinetic energy spectra of flaton (left) and AD (right) sectors, thermally nor- 
malized (top) so that a thermal spectrum is flat, and energy normalized (bottom) so that 
the area under the graph is the total energy (CP— ): |0a;P/Vo, \'ipk\'^ /V2, /c^|0a:P/'"^pqVo, 
k^\i'k?/mlHy2 at t = 0, 40, 80, , 160, 200, 250, 300, 500, 1000 versus k/m. Relevant 
mass scales are mpg = 0.11m, m^ = 0.47m, m^H^ = 0.61m, m^ ^ = 0.73m, 



mjT u = 1.8 m. 



very efficient. As can be seen from the top right graph, all the AD lattice modes have 
thermalized by about t ~ 250 m~^. Thus the AD preheating is artificially truncated by 
our lattice cutoff, so in reality should be even more efficient, bringing the AD fields in 
even closer to the origin. However, the preheating available with our limited lattice is 
sufficient to lead to excellent conservation of lepton number, as shown in Figure 12. 

The energy densities in the flaton sector are shown in Figure 17. The flrst flaton 
oscillation is dominated by radial kinetic and potential energy, as would be expected for 
an essentially homogeneous oscillation. During the second oscillation, a large part of 
the potential energy is converted into angular gradient and then angular kinetic energy, 
causing the radial flaton to settle down around its minimum. The late time axion 
dispersion oscillations are again evident from the angular kinetic and gradient energies. 

The energy densities in the AD sector are shown in Figure 18. Initially / is sitting in 
its thermal inflation minimum and hd is at the origin, hence the potential energy is —Vi 
and the kinetic and gradient energies are negligible. Then the AD potential gets lifted 
up by the flaton and the AD flelds roll in towards, and oscillate homogeneously about, 
the origin with negligible gradient energy. As the flaton returns near the origin, the 
AD potential is dropped and / returns to its thermal inflation minimum, though now 
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Figure 17: Kinetic, gradient and potential energy densities of the radial (solid) and 
angular (dotted) flaton fields averaged over the lattice as a function of time (CP—): 
P(j)/Vo versus mt. 
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Figure 18: Kinetic, gradient and potential energy densities of the AD sector averaged 
over the lattice as a function of time (CP—): Pad/V2 versus mt. 
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Figure 19: Radial (solid) and angular (dotted) energy densities of the fiaton averaged 
over the lattice as a function of time {CP—): p^JVq versus mt. 
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Figure 20: Fiaton and AD energy densities averaged over the lattice as a function of 
time [CP—]: P(j,/Vq, Pad/Vq versus mt. 
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oscillating with substantial kinetic energy. The flaton then rolls out again and settles 
around its minimum, permanently lifting the AD potential. The AD fields once again 
roll in towards, and oscillate about, the origin, though this time with rapidly growing 
gradient energy. Finally, the AD energies all gradually increase, though at a decreasing 
rate, due to energy transfer from the fiaton sector. 

Figure 19 shows the energy transfer between the radial and angular fiaton sectors. 
The fiaton energy is initially all in the radial component but quickly becomes equally dis- 
tributed between the radial and angular components. Over much longer timescales, far 
beyond the scope of this simulation, the Hubble expansion, and possibly also preheating 
and decay, will reduce the amplitude of the flaton oscillations about its minimum. Then 
we expect the radial fiaton and the axion to decouple and the energy in the axion sector 
to red-shift away, leaving the universe dominated by the radial fiaton until it finally 
decays. However, the details of this process deserve further investigation. 

Figure 20 shows the energy transfer from the fiaton sector to the AD sector. Since 
the fiaton sector dominates, Vq ^ V^i, V2, the energy transfer is negligible for the fiaton 
but is significant for the AD sector. The energy transfer could be dangerous for the 
conservation of the lepton number if it increases the amplitude of the / field too much. 
Whether this happens depends on the efficiency of the energy transfer from the fiaton 
sector to the AD sector, which from Figure 20 seems to be rapidly declining, the energy 
fiow within the AD sector to higher k modes, which from Figure 16 seems to be very 
efficient but is artificially cut off at an early stage in our simulation due to the finite 
lattice size, and the transfer of energy from the AD sector to other sectors via thermal 
friction [b], which is beyond the scope of this paper. However, as the lepton is well 
conserved in our simulation, see Figure 12, despite artificially cutting off the preheating 
and not including thermal friction at all, we are confident it will be conserved in reality. 

Thus, despite the dynamics being complex, our simulation suggests that our baryo- 
genesis mechanism does work in this model. 

5 Conclusion 

In this paper, we proposed a cosmological model based on a simple extension of the 
MSSM with superpotential 



W = XuQH^u + XdQHj+ X.LHae + X^^'^H^H^ + -X, {LH.,f + X^<t>xx (240) 



and key parameter condition Eq. (28). The model is obtained from our original thermal 
infiation and baryogenesis model of Ref. [7] by removing the fiaton self-interaction term 
A^c/)^, with the fiaton now stabilized by the renormalization group running of its poten- 
tial. Removing the fiaton self-interaction gives the model a PQ symmetry and axion, 
with the scales required for thermal infiation and axions well matched. As shown in 
Section 4, our baryogenesis scenario works well in this model, as it did in our original 
model [8], but now in an even more minimal setting, using essentially all the terms in 
the potential. 
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Comparison with multi-field thermal infiation axion models 

Our model, though simpler, has much in common with the multi-field thermal inflation 
axion models studied in Refs. [48, 49, 50, 51], and our baryogenesis scenario would 
probably work in those models too, but there are at least two important differences. 
First, in the absence of the flaton self-interaction, the flatino/axino becomes light making 
it the LSP, see Eq. (56). Second, the renormalization group running stabilization of 
the flaton's potential suppresses the flaton's mass squared around the minimum of its 
potential, see Eq. (49). 

An axino LSP would be over-produced in the standard hot Big Bang cosmology, but 
in our model the flaton decays late at a temperature of order 1 GeV. The lack of a 
flaton self-interaction suppresses direct axino production in the flaton decay, but axinos 
are still typically over-produced either in the flaton decay or by the decay of NLSPs 
in the thermal bath. This provides the strongest constraint on our model, forcing the 
axino mass to be at the lower end of its expected range, rria ~ 1 GeV. Avoiding over- 
production of axinos also requires /„ > 10^^ GeV, which is nonetheless nicely consistent 
with the scale required for thermal inflation, our baryogenesis scenario and axion dark 
matter, see Section 3. 

The suppressed flaton mass suppresses the flaton decay rate to hot axions, whose 
over-production would conflict with Big Bang nucleosynthesis, and is the strongest con- 
straint on the multi-fleld thermal inflation axion models [48, 49, 50, 51]. It also tends 
to put the flaton mass below threshold for Higgs production, reducing the decay tem- 
perature and helping avoid over-production of dark matter. 

Observational tests and signatures 

Our model has a variety of observational tests and signatures. Thermal inflation wipes 
out any gravitational waves at solar system scales generated during primordial inflation 
[62], so observation of such would rule out thermal inflation. However, the flrst order 
phase transition that ends thermal inflation generates its own gravitational waves [57]. 
These may be observable at future gravitational wave detectors such as BBO and DE- 
CIGO, and a correlated analysis of ultimate-DECIGO [■ ] may probe the heart of the 
thermal inflation parameter space. 

Thermal inflation and the following period of flaton domination redshift primordial 
perturbations by 10 to 15 e-folds compared with a standard hot Big Bang history. This 
may have an observable effect on the primordial density perturbations. For example, 
many simple models of primordial inflation predict a spectral index n — 1 = —a/N + 
0{1/N^), where N is the number of e-folds between horizon exit and the end of the 
primordial inflation. Thermal inflation would reduce N by 10 to 15. One could also try 
to reconstruct A^. For example, in the simple class of models above, (n — l)^/n' = —a 
determines the model, {n — l)/n' = N determines the number of e-folds, and {n — 
l)n" /{n'Y = 2 provides a check of assumptions, where n' = dn/d\nk, etc. 

The key parameter condition for our baryogenesis scenario, Eq. (28), can be tested 
at future accelerators. The roll away of the AD fleld may extend the thermal inflation 
by 5 or 6 e-folds, see Eqs. (122) and (130). The x and x flelds, needed to couple the 
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flaton to the thermal bath, acquire intermediate scale masses after thermal inflation 
and, assuming that they are not all MSSM singlets, will affect the renormalisation of 
the MSSM couplings. The late decay of the fiaton may dilute, or even enhance, the dark 
matter abundance compared with a standard hot Big Bang history. 

Thermal inflation axion models can be tested by the effect of the hot axions on Big 
Bang nucleosynthesis. Significant parts of the parameter space of multi-field thermal 
inflation axion models are already ruled out by this test [18, 49, 50, 51]. As discussed in 
Sections 2.5.1 and 3.7, and above, the hot axion production is suppressed in our model 
but is still expected to be greater than the usual thermal production, so this might also 
provide a test of our model in the future. 

Our model, and Moxhay and Yamamoto's original single-field flaton axion model 
[55], can be tested by its prediction that the LSP is the axino/flatino. Furthermore, in 
our model we expect rria ~ 1 GeV, see Sections 2.2 and 3.4. This provides an important 
test of our model at future accelerators [64, 65, 66, 67]. Also, the tightness of the axino 
bounds, see Figures 9 and 10, suggests that the dark matter should be composed of 
significant amounts of both axions and axinos. 

Summary 

In summary, our simple and natural extension of the MSSM leads to a rich but remark- 
ably consistent cosmology combining thermal inflation, baryogenesis, axions and axinos, 
with observational implications for primordial perturbations, gravitational waves. Big 
Bang nucleosynthesis, dark matter and particle accelerators. 
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A Flaton decay rates 



A.l Decay to axions 

Decomposing the flaton as in Eq. (2), the flaton kinetic term generates the radial flaton 

axion interaction 

5r {day 



v^0o 



Since rua ~ 0, we get the flaton decay rate to axions 

,3 



r. 



m 



PQ 



647r0^ 



(241) 



(242) 



A. 2 Decay via SM Higgs 

The superpotential term X^cp'^HuHa generates the radial flaton SM Higgs interaction 



iA,n0i^ 



Hu\' + W) -. V2 



l/^l 



6r h'^ 
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where 



I R|2|,,|2 



^\ = -"^k+"^L + 2|/^r 



I R|2\ I ,,|2 



m\ ) 0, 



(243) 
(244) 



and we have used Eq. (103). Therefore the flaton decay rate via (possibly virtual) SM 
Higgses is [68] 
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(246) 
iPQ'/^o V "^~aJ 

and F/j is the decay rate of the SM Higgs. The dominant decay mode of the SM Higgs 
is to bottom quarks, giving [21] 
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(247) 



where g2 — 0.65 is the SU(2) gauge coupling, m;, ~ 5 GeV is the bottom quark mass 
and my/ ~ 80 GeV is the W boson mass. Therefore, defining x = p/mpq, y = q/mpQ 
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and $, = rrih/mpQ, we have 
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(249) 



A. 3 Decay to h and h via flaton-Higgs mixing 

The superpotential term X^cjP'HuHd, induces flaton-Higgs mixing, allowing the flaton 
mass eigenstate to decay directly to SM fermions. 
The neutral CP even Higgs states are 



Hi = el.oh + e^oH 

iJ° = Ql,oh + Q^oH 



(250) 
(251) 



where h and H are the mass eigenstates, h being the SM Higgs. The mass matrix 
elements are 



Ml 



-m\ (1 + cos 2/3) + -m| (1 - cos 2/3) 



M^ 






MjjOjjO 



{m\ + m|) sin 2/3 



-m\ (1 - cos 2/3) + -m| (1 + cos 2/3) 



giving the tree level Higgs masses 
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where 



sin 2/3 



2 I 2 I o I 1 2 

2^ 
mi 

' tan2/3-l ' ^'^' 

47 



(252) 
(253) 
(254) 

(255) 
(256) 

(257) 
(258) 

(259) 



A is the neutral CP odd Higgs boson and tan/3 is defined as 

tan/3=— (260) 

Vd 

where Vu and Vd are the magnitudes of the vacuum expectation values of Hu and Hd, 
and the electroweak symmetry breaking scale is 

v^=vl + vj = {174GeYf (261) 

The flaton mixes with h and H 

h = Qlh + Qf^H + ei^ (262) 

H = Q%h + eiH + 0^4) (263) 

= ej/i + ef^ + ej0 (264) 

where h, H and are the mass eigenstates, due to the mass matrix elements 

Ml^ = m^Q (265) 

Mlo^ = 2 sin f3 (2\fi\^- ml cos^/3)^ (266) 

Mlo. = 2cos/?(2|/i|2-m^sin2/5) ^ (267) 

Since the decay is dominated by the bottom quark channel for m^ < 1 TeV [ ] 

XbH^bb -^ Xb&io^bb (268) 
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we are only interested in the H^ mixing 

<o = e^0f + e^oei (270) 
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Diagonalizing the Higgs mass matrix gives 
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(271) 

(272) 
(273) 

(274) 



and perturbatively diagonalizing the flaton Higgs mass matrix gives 
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et 



e 



H 



^loHoMlo, 



A//2 



m. 



^m 



(m 



Ml,)MMl^.^.^-ml) +M 






4 1- 



\Bl 



lyUpf 



{ml - ml^) 00 



A© + 5, 



e 



4|/i|2-2m^cos2/3 



sin/3 



(^is 



H'l 



mi 



v/A| + 2Ae5esin/5 + S| 



("^1. - -^^J V (-^Ho^o - m|) + A^^^o 



m\v sin 4/? 



2 (m2^ - m2 q) 



Ae + 5, 



e 



2 
A 



v/A| + 2Ae5esin/5 + 5, 



Be 



{m.\ + m|) sin/5 
m| cos 2/5 — m^ 



(275) 



(276) 



(277) 



(278) 



(279) 
(280) 



Thus the mixing coefficient 0*^0 has a comphcated dependence on the parameters, but 
can be roughly parameterized as 
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with Cq ~ 1. With this parameterization, Eq. (269) becomes 
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where A is defined in Eq. (246). 



A. 4 Decay to gluons via x ^^<^ X 

The superpotential term \(pxx allows the fiaton to decay to SM gauge fields via a 
virtual XX loop. The dominant decay channel is via heavy quark superfields to gluons 



«i^'"^PQ f^ ^ 95«3 



'^^ 1447r3(/.2 V" ' 47r 



(283) 



where as = gl/{A7r) ~ 0.1. 
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A. 5 Decay to axinos 

Axino mass is generated at one-loop and given by Eq. (56) 



^^ = iE4^x (284) 



167r2 

X 

The 101 dependent renormalisation of A^^, and A^ is [69] 

5^ = |5(ElVP + 2|A/-45:a(x)«J) (285) 

r^ A 

^ • ' ■ -■■ - ■ • — " 286) 



din 
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where the Ca{x) are the quadratic Casimir invariants, for example C^ix) = 4/3 for a 
quark superfield x- This generates the coupling of the flaton to the axino, given in 
Eq. (136), with 

dlnrria , . 

Sx^xfex'IVl'^x' + E,'|Axf^x + 4|A,rA, + 4E.a(x)^^(M„-A, 



(288) 
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